Monise Terra Cerezini

Exercicios da aula 4

#Cervejas

Uma amostra de 30 estudantes foi indagada sobre seu tipo de cerveja preferida, com o seguinte resultado:

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")

1. Represente este resultado como um grfico de barras e um dotplot (funo dotchart). Qual tem maior razo dado/tinta? 

barplot(table(cervejas))
dotchart(table(cervejas))

#Rios

O conjunto de dado rivers est em um objeto do pacote datasets. e tem a extenso dos principais rios da Amrica do Norte, em milhas (consulte a ajuda para detalhes). Para carregar um objeto de dados de um pacote em sua rea de trabalho use a funo data: 

data(rivers)

1. Qual a proporo do total de rios que tem extenso menor do que a mdia?

>mean(rivers<mean(rivers))
[1] 0.6666667

2. Qual  o quantil de 75%?

>quantile(rivers, prob=.75)
75% 
680 

3. Compare a mdia, mdia truncada a 25% e mediana. H diferenas?

> mean(rivers,trim=25)
[1] 425
> mean(rivers)
[1] 591.1844
> median(rivers)
[1] 425


#Caixetais

caixeta <- read.csv("caixeta.csv", as.is=T)
str(caixeta)

1. Construa um histograma da altura dos fustes do caixetal.

hist(caixeta$h)

2. Construa histogramas da altura das rvores para os diferentes caixetais ('local').

hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="retiro"])
hist(caixeta$h[caixeta$local=="jureia"])

3. H diferenas entre as estruturas (distribuio de tamanhos) dos caixetais? 

par(mfrow=c(1,1))

histogram(~caixeta$h|caixeta$local) #No parece haver diferena entre os caixetais


#Eucaliptos

1. Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regio ('regiao') e rotao ('rot'). 

eucaliptos <- read.csv("egrandis.csv", as.is=T)
str(eucaliptos)

library(lattice)
bwplot(dap~regiao|rot,data=eucaliptos,ylab = "Dap")

2. Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.

qqnorm(eucaliptos$ht)
qqline(eucaliptos$ht)


#Mais Caixetais

1. Analise a relao dap-altura ('dap' e 'h') em funo do caixetal, mas somente para os fustes de caixeta (Tabebuia cassinoides). 

caixeta.tabebuia=caixeta[caixeta$especie=="Tabebuia cassinoides",]
str(caixeta.tabebuia)

plot(caixeta.tabebuia$cap,caixeta.tabebuia$h,ylab="DAP", xlab="h")

2. Para a mesma relao do item anterior, verifique linearidade com a funo scatter.smooth

scatter.smooth(caixeta.tabebuia$cap~caixeta.tabebuia$h,ylab="DAP", xlab="h")

3. Utilizando o pacote lattice, analise a relao dap-altura ('dap' e 'h') em funo do caixetal, mas somente para os fustes de caixeta (Tabebuia cassinoides).

xyplot(cap~h|local,data=caixeta.tabebuia,ylab="DAP", xlab="h")
